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Abstract 

Assuming the Majorana nature of massive neutrinos, we generalize the 
Friedberg-Lee neutrino mass model to include CP violation in the neutrino 
mass matrix M v . We show that a favorable neutrino mixing pattern (with 
9 12 ps 35.3°, 6» 23 = 45°, 13 / 0° and <5 = 90°) can naturally be derived from 
M v , if it has an approximate or softly-broken fj,-r symmetry. We point out a 
different way to obtain the nearly tri-bimaximal neutrino mixing with 5 = 0° 
and non- vanishing Majorana phases. The most general case, in which all the 
free parameters of M v are complex and the resultant neutrino mixing matrix 
contains both Dirac and Majorana phases of CP violation, is also discussed. 
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I. INTRODUCTION 



The solar [1], atmospheric [2], reactor [3] and accelerator [4] neutrino experiments have 
provided us with very convincing evidence that neutrinos are massive and lepton flavors 
are mixed. Given the basis in which the flavor eigenstates of charged leptons are identified 
with their mass eigenstates, the phenomenon of neutrino mixing can simply be described 
by a 3 x 3 unitary matrix V which transforms the neutrino mass eigenstates u 2 , v 3 ) into 
the neutrino flavor eigenstates (u e , v , v T ). We assume massive neutrinos to be Majorana 
particles and parametrize V as 
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sin 9ij (for ij = 12,23 and 13). A global analysis of current 
experimental data [5] yields 30° < 9 12 < 38°, 36° < 9 23 < 54° and 9 13 < 10° at the 99% 
confidence level, but three CP- violating phases of V (i.e., the Dirac phase S and the Majorana 
phases p and a) are entirely unrestricted. In order to interpret the largeness of two neutrino 
mixing angles together with the smallness of three neutrino masses, many theoretical and 
phenomenological models of lepton mass matrices have been proposed in the literature [6] . 
Among them, the scenarios based on possible flavor symmetries are particularly simple, 
suggestive and predictive. 

In this paper, we focus our interest on the neutrino mass model proposed recently by 
Friedberg and Lee (FL) [7]. The neutrino mass operator in the FL model is simply given by 

^-mass = a(y T - 17 J { V r ~ ^) + b ~ V e) ~ ^e) + C fa - V T ) (u e - U T ) 

+ m o(¥e + Vf +I7 A) , (2) 

where a, b, c and m are all assumed to be real. A salient feature of C u _ mass is its partial 
gauge-like symmetry; i.e., its a, b and c terms are invariant under the transformation v a — > 
v a + z (for a = e, p, r) with z being a space-time independent constant element of the 
Grassmann algebra [7]. Eq. (2) means that the neutrino mass matrix M v takes the form 
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Diagonalizing M v by the transformation ^pl^^fl = Diag{m 1 , m 2 , m 3 }, in which ra i (for 
i = 1, 2, 3) stand for the neutrino masses, one may obtain the neutrino mixing matrix 
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where 9 is given by tan 9 
following observations: 
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c) / [(b + c) — 2a] . This interesting result leads us to the 
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• If 9 = 0° holds, V FL will reproduce the exact tri-bimaximal neutrino mixing pattern 
(with tan# 12 = l/y/2 or 12 ss 35.3°, 23 = 45° and 13 = 0°) [8]. The latter, which can 
be understood as a geometric representation of the neutrino mixing matrix [9], is in 
good agreement with current experimental data. Non-vanishing but small 9 predicts 
sin0 13 = (2/v / 6)sin(0/2), implying 9 < 24.6° for 13 < 10°. On the other hand, 9 23 
will mildly deviate from its best-fit value 6 23 = 45° if 9 (or 13 ) takes non-zero values. 

• The limit = 0° results from b = c. When b = c holds, it is easy to check that the 
neutrino mass operator £„_ mass has the exact p-r symmetry (i.e., £„_ mass is invariant 
under the exchange of p and r indices). In other words, the tri-bimaximal neutrino 
mixing is a natural consequence of the p-r symmetry of M v in the FL model. Then 
9 13 7^ 0° and 9 23 ^ 45° measure the strength of p-r symmetry breaking, as many 
authors have discussed in other neutrino mass models [10]. 

In addition, one may consider to remove one degree of freedom from C u _ mass or M v (for 
instance, by setting c = [7]). 

We aim to generalize the FL model to accommodate CP and T violation for massive 
Majorana neutrinos L The effective Majorana neutrino mass term can be written as 



£™ = 2 ("» "„> L M »\»lj + h - C - > ( 5 ) 

where v c a = Cz/J (for a — e, p, r), and M v is of the same form as that given in Eq. (3). Now 
the parameters of M v (i.e., a, b, c and m ) are all complex. Then we are able to derive both 
the Dirac phase (S) and the Majorana phases (p and a) for the neutrino mixing matrix V. 
Two special cases are particularly interesting: 

• Scenario (A): a and m Q are real, and b = c* are complex. We find that the /i-t 
symmetry of M v is softly broken in this case, leading to the elegant predictions 9 13 ^ 0°, 
#23 = 45° and 8 = 90°. Two Majorana phases p and a keep vanishing. 

• Scenario (B): a, b and c are real, but m is complex. We find that the results of 9 12 , 
9 23 and 9 13 obtained from V FL keep unchanged in this case, but some nontrivial values 
of the Majorana phases p and a can now be generated. The Dirac phase 5 remains 
vanishing. 

The remaining part of this paper is organized as follows. More detailed discussions 
about scenarios (A) and (B) will be presented in section II. Section III is devoted to a 
generic analysis of the generalized FL model with no special assumptions. Finally, we make 
some concluding remarks in section IV. 



l To include T violation into the model, Friedberg and Lee [7] have inserted the phase factors e ±i?? 
into Eq. (2) by replacing the term c (y e — V T ) {y e — u T ) with the term c {e~ %rt V e — V T ) {e rlT] v e — u T ). 
The resultant neutrino mass matrix is no longer symmetric, hence it definitely describes Dirac 
neutrinos instead of Majorana neutrinos. 
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II. TWO SIMPLE SCENARIOS 



First of all, let us consider two special but interesting scenarios of the generalized FL 
model and explore their respective consequences on three neutrino mixing angles and three 
CP-violating phases. They are quite instructive in phenomenology and may easily be tested 
by a variety of long-baseline neutrino oscillation experiments in the near future. 



A. Scenario (A) 



In this scenario, a and m are real, b and c are complex, and b = c* holds. The corre- 
sponding neutrino mass matrix M v reads 
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We find that M v can be diagonalized by the transformation V^MJV* = Diag{m 1; m 2 , m 3 }, 
where 
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and S is just the Dirac phase of CP violation defined in the standard parametrization of V. 
A straightforward calculation yields 5 = 90°, 



tan# = 



V3 Im 0) 



m + a + 2Re(6) ' 
together with three neutrino masses 

m 1 = ^/[m + a + 2Re (b)] 2 + 3 [Im (b)] 2 - a + Re (b) , 
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Note that the special result of 5 is a natural consequence of the purely imaginary term 
b — c = 2ilm (b) in this scenario. The difference between b and c can be referred to as the 
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We observe that the only difference between V in Eq. (fO) and V FL in Eq. (4) is the 
introduction of a special CP- violating phase (see also Ref. [11] for a discussion about the 
maximal leptonic CP violation with 8 = 90°). This CP- violating phase, which is attributed 
to the soft breaking of /i-r symmetry, can change the prediction of V FL for 9 23 . Comparing 
Eq. (10) with Eq. (1), we immediately obtain 2 

S in #12 = 7k , X Q , 
V 2 + cos 9 

1 

72' 
2 

sin #i3 = ^ sin 2 ' ( n ) 

together with <5 = 90° and p = a = 0°. Without loss of generality, we have restricted 6 to 
the first quadrant. The leptonic Jarlskog parameter J [12], which is a rephasing-invariant 
measure of CP violation in neutrino oscillations, reads J = sin#/(6y / 3)- One can see that 
the soft breaking of /i-r symmetry leads to both 9 13 ^ 0° and J ^ 0, but it does not affect 
the favorable result 9 23 = 45° given by the tri-bimaximal neutrino mixing pattern. On the 
other hand, sin 9 12 ~ 1/ y/3 is an excellent approximation, since 9 must be small to maintain 
the smallness of 9 13 . In view of 9 13 < 10°, we obtain 9 < 24.6° and J < 0.04. It is possible 
to measure J ~ (9(10~ 2 ) in the future long-baseline neutrino oscillation experiments. 

Note that the neutrino masses rely on four real model parameters m , a, Re (6) and Im(6). 
Thus it is easy to fit the neutrino mass-squared differences Am^ = (7.2 • • -8.9) x 10~ 5 eV 2 
and Am 2 2 = ±(2.1 •••3.1) x 10~ 3 eV 2 [5]. Such a fit should not involve any fine-tuning, 
because (a) the number of free parameters is larger than the number of constraint conditions 
and (b) three neutrino masses have very weak correlation with three mixing angles. A 
detailed numerical analysis shows that only the normal mass hierarchy (m 1 < m 2 < m 3 ) is 
allowed in this scenario. FIG. 1 illustrates the parameter space of m , a, Re (6) and lm(b), 
where m < 0.2 eV has typically been taken as a generous upper bound on the absolute 



neutrino mass [13]. Because of m = m 2 , it is straightforward to get m > y Am 21 m 
0.009 eV, as shown in FIG. 1. The small fx-r symmetry breaking (or small 9) requires 
the small magnitude of Im(fo). Thus Eq. (10) allows us to get Re (6) (m 1 — m 2 )/3 
in the Im(6) — > limit, implying that Re (6) is negative and its magnitude is very small. 
Furthermore, a ~ (m 3 — m 2 )/2 approximately holds in the Im(6) — ■> limit, implying that 



the maximal value of a is roughly half of m 3 ~ Y|Am 2 2 | (i.e., a < 0.025 eV) when three 
neutrino masses develop a strong normal hierarchy. 

B. Scenario (B) 

In this simple scenario, only m is assumed to be complex. While the neutrino mass 
matrix M v takes the same form as that given in Eq. (3), its diagonalization (V* M U V* = 



2 It is trivial to redefine the phases of charged-lepton fields, such that the location of CP-violating 
phases in Eq. (10) is the same as that in Eq. (1). 
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where (f) i (for i = 1,2,3) originate from the imaginary part of to . Comparing between 
Eqs. (1) and (12), one can see that the Majorana phases of CP violation in the standard 
parametrization of V are given by p = <p 1 — 3 and a = <f) 2 — <f) 3 (i.e., the phase factor e*^ 
in Eq. (12) is finally rotated away by redefining the phases of three charged- lepton fields). 
After a straightforward calculation, we obtain 
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The final result for the neutrino mixing matrix V is 
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where p and a are given by 
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One can see that the only difference between V in Eq. (16) and V Fh in Eq. (4) is the 
introduction of two Majorana phases of CP violation. Although p and a have nothing to 
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do with the behaviors of neutrino oscillations, they may significantly affect the neutrinoless 
double-beta decay [14]. Comparing Eq. (16) with Eq. (1), we arrive at 

1 
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yjl + cos 9 
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together with 5 = 0° for the Dirac phase of CP violation. Again 9 has been restricted to the 
first quadrant. The results for 9 12 and 9 13 in this scenario are the same as those obtained in 
scenario (A), but the Jarlskog parameter J is now vanishing. Because of the /i-t symmetry 
breaking, 9 23 may somehow deviate from the favorable value 9 23 = 45°. Given 9 < 24.6° 
corresponding to 9 13 < 10°, 9 23 is allowed to vary in the range 37.8° ^ 9 23 < 45° [15]. 



Note that the neutrino masses depend on five real model parameters a, b, c, Re( 



and Im(m ). Hence there is sufficient freedom to fit two observed neutrino mass-squared 
differences Am^ and Am| 2 . A detailed numerical analysis shows that both the normal 
mass hierarchy (m 1 < m 2 < m 3 ) and the inverted mass hierarchy (m 3 < m 1 < m 2 ) are 
allowed in scenario (B). FIGs. 2 and 3 illustrate the parameter space of a, b, c, Re(m ) 
and Im(m ), where |m | < 0.2 eV has been taken [13]. Taking account of the analytical 
reciprocity between b and c (i.e., an exchange of b and c in Eq. (15) does not affect the 
result of m i obtained in Eq. (14)), we plot their allowed regions in the same figure without 
any confusion. Note that m cannot be purely imaginary in both normal and inverted 
mass hierarchies, otherwise we would be left with Am 2l < 0, which is in disagreement with 
current experimental data. Moreover, the lower bound of |Re(m )| is obviously restricted 



by the magnitude of m 2 (i.e., m 2 ~ J Am 2l m 0.009 eV in the normal mass hierarchy or 



m 2 ~ ^\Am 32 \ ~ 0.05 eV in the inverted mass hierarchy). In view of (m + — m_) oc {m\—m\) 
from Eq. (14), we find m+ ^ m_. This inequality implies that a = b = c is not allowed, 
as one can see from Eq. (15) or FIGs. 2 and 3. The disconnected regions in the plots 
of a versus b (or c) are ascribed to the ambiguity induced by the sign of Re(m ). FIG. 4 
shows the allowed regions of p and a. We see that both of them are less restricted, as a 
consequence of the large freedom associated with the imaginary part of m . 



III. GENERIC ANALYSIS 

Now let us assume all the parameters of M v in Eq. (3) to be complex and calculate its 
mass eigenvalues and flavor mixing parameters. Using the tri-bimaximal mixing matrix [8] 
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we transform M v into the following form: 
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This complex mass matrix can be diagonalized by the transformation V±(VqM v Vq)V{ = 
Diag{m 1 , m 2 , m 3 }, where 
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The neutrino mixing matrix V turns out to be V = V V 1 , in which 5 is just the Dirac phase 



of CP violation. Of course, p 
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are the Majorana phases of CP 



violation in the standard parametrization of V. 

To be explicit, we express 9 and 5 in terms of the parameters of M u . The results are 
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Furthermore, three mass eigenvalues of M v and two Majorana phases of V are found to be 
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The final result of V is 
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from which we obtain the Jarlskog parameter J = sin#sin<5/(6\/3)- Comparing Eq. (27) 
with Eq. (1), we arrive at 
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One can see that the predictions for 9 12 and 9 l3 are the same as those obtained in scenarios 
(A) and (B). Thus they are two typical and general consequences of the FL model. As for 
#23, the generic result in Eq. (28) may easily reproduce the special result in Eq. (11) for 
scenario (A) with 5 = 90° or that in Eq. (18) for scenario (B) with 5 = 0°. Given 9 < 24.6° 
(i.e., 9 13 < 10°), 6> 23 is allowed to vary in the range 37.8° < 9 23 < 52.2° for arbitrary 5. 

As m , a, b and c are all complex, we now have much more freedom to fit two observed 
neutrino mass-squared differences Am^ and Am 32 . Of course, both the normal neutrino 
mass hierarchy (m 1 < m 2 < m 3 ) and the inverted neutrino mass hierarchy (m 3 < m 1 < m 2 ) 
are expected in this generic case. We shall not carry out a numerical analysis of the parameter 
space, however, just because it involves a lot of uncertainties and is not as suggestive as that 
in scenario (A) or scenario (B). 



IV. CONCLUDING REMARKS 

We have pointed out that the nearly tri-bimaximal neutrino mixing can be regarded 
as a natural consequence of the slight /i-t symmetry breaking in the FL neutrino mass 
model. Another straightforward consequence of the /i-t symmetry breaking is CP violation. 
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Assuming the Majorana nature of massive neutrinos, we have generalized the FL model to 
introduce the CP-violating effects. In addition to a generic analysis of the generalized FL 
model, two simple but intriguing scenarios have been proposed: scenario (A) involves the 
softly-broken p-r symmetry, leading to the elegant predictions 9 l3 ^ 0°, 9 23 = 45°, 5 = 90° 
and vanishing Majorana phases of CP violation; scenario (B) predicts 9 13 ^ 0°, 9 23 ^ 45°, 
5 = 0° and two nontrivial Majorana phases of CP violation, on the other hand. Both 
scenarios are compatible with current experimental data. 

Although our discussions about the generalized FL model are restricted to low-energy 
scales, it can certainly be extended to a superhigh-energy scale (e.g., the GUT scale or the 
seesaw scale). In this case, one should take into account the radiative corrections to both 
neutrino masses and flavor mixing parameters when they run from the high scale to the 
electroweak scale. A particularly interesting point is that the Majorana phases p and a can 
radiatively be generated from the Dirac phase 5 in scenario (A), while the Dirac phase 5 can 
radiatively be generated from the Majorana phases p and a in scenario (B) [16]. Thus the 
features of leptonic CP violation can fully show up in both scenarios at low-energy scales, if 
they are originally prescribed at a superhigh-energy scale (see Ref. [17] for a detailed analysis 
of the running behaviors of CP- violating phases in the nearly tri-bimaximal neutrino mixing 
pattern, including the minimal supergravity threshold effects). 

We conclude that the p-r symmetry and its slight breaking are useful and suggestive for 
model building. We expect that a stringent test of the generalized FL model, in particular 
its two simple and instructive scenarios, can be achieved in the near future from the reac- 
tor neutrino oscillation experiments (towards measuring 9 l3 ), the accelerator long-baseline 
neutrino oscillation experiments (towards detecting both 9 l3 and 5) and the neutrinoless 
double-beta decay experiments (towards probing the Majorana nature of massive neutrinos 
and constraining their Majorana phases of CP violation). 
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FIG. 1. The parameter space of (m , a) and (Re(6), Im(fe)) in Scenario (A), where only the 
normal neutrino mass hierarchy is allowed. 
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FIG. 2. The parameter space of (Re(m ), Im(m )) and (a, 6 or c) in Scenario (B) with the 
normal neutrino mass hierarchy. 
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FIG. 3. The parameter space of (Re(m ), Im(m )) and (a, 6 or c) in Scenario (B) with the 
inverted neutrino mass hierarchy. 
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